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We discuss the general transport properties of superconducting quantum point contacts. We 
show how these properties can be obtained from a microscopic model using nonequilibrium Green 
function techniques. For the case of a one-channel contact we analyze the response under different 
biasing conditions: constant applied voltage, current bias and microwave-induced transport. Current 
on : fluctuations are also analyzed with particular emphasis on thermal and shot-noise. Finally, the case 

of superconducting transport through a resonant level is discussed. The calculated properties show a 
remarkable agreement with the available experimental data from atomic-size contacts measurements. 
We suggest the possibility of extending this comparison to several other predictions of the theory. 

'd 

^ ■ I. INTRODUCTION 

<D . 

Since the discovery of the Josephson effect [|lj the electronic transport between weakly coupled superconducting 
electrodes (weak superconductivity) has been a subject of growing interest Q. Typically, weak superconductivity has 
been studied in SIS, SNS and S-c-S junctions, where S, N, I and c denote superconductor, normal metal, insulator 
and constriction respectively. Recent technological advances have made possible the fabrication of mesoscopic S-c-S 
junctions in which the electrodes are connected by a small number of conduction channels. These systems are usually 
referred to as superconducting quantum point contacts (SQPC), examples of which are the S-2DEG-S junctions ||] 
and atomic contacts produced by break junctions [|]J^] and scanning tunneling microscope (STM) || techniques. 

On the theoretical side there has also been a parallel advance with the development of fully quantum mechanical 
theories for the transport properties of superconducting one-channel contacts Jz| |l0|] • There has been a remarkable 
y—i , agreement between theoretical predictions and experimental results for the quantities that have so far been measured. 
These quantities include the phase-dependent supercurrent in a high transmissive contact p"l| ] and the dc current at 
constant bias voltage |^,^]. As we discuss in this paper, there remain many exciting predictions of the microscopic 
' theories to be explored experimentally. 
_ , The aim of this paper is to present an overview of the main theoretical results that have been obtained for different 
t-H ' microscopic models of an SQPC. An interesting aspect of superconducting transport is that qualitatively different 
behaviors are exhibited depending on how the system is biased. This will be analyzed in this work by discussing 
the cases of phase, voltage and current bias together with the case of transport under microwave radiation. The 
models are introduced in section II together with the nonequilibrium Green functions formalism used to calculate 
their transport properties. Section III is devoted to the voltage biased case for which we discuss the comparison of the 
fully quantum mechanical calculation with semiclassical standard theories and the available experimental results. We 
also discuss the limit of very small voltage. In section IV the current biased case is briefly analyzed while the response 
under microwave radiation and its possible relevance for directly detecting Andreev states is discussed in section V. 
Thermal and shot-noise are the subject of section VI where we discuss the conditions for observing coherent transport 
of multiple charge quanta from the noise-current ratio. Finally in section VII the superconducting transport through 
a resonant level is analyzed both in the limits of very large and very small charging energy. The general conclusions 
are summarized in section VIII. 



II. MICROSCOPIC MODEL AND GREEN FUNCTION FORMALISM 



A schematical representation of a quantum point contact is depicted in Fig. 1 . For a typical point contact the length 
of the constriction between the electrodes, Lc, is much smaller than the superconducting coherence length £o and 
its width Wc is ~ Xf, the electron Fermi wavelength. The first condition ensures that the detailed superconducting 
phase and electrochemical potential profiles in the constriction region are unimportant and can be safely approximated 
by step functions. On the other hand, the condition W ~ Xf implies that there are only a few conduction channels 
between the electrodes. 
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FIG. 1. Schematic representation of a superconducting quantum point contact. 



The general mean-field Hamiltonian for a superconducting system can be written in terms of the electron field 
operators ipa{r) 



6 = J df 



r) + A* {r)i,\ {r)i>\ (r) + A(r> x (r)V T (r) j (1) 

where 7i e is the one-electron Hamiltonian and A(r) is the superconducting order parameter. The problem of calcu- 
lating transport properties in such a continuous representation for a non-homogenuous system is extremely involved 
requiring the knowledge of the adequate boundary conditions at the interfaces. Some attempts in this direction have 
been recently presented by Zaitsev and Averin []12| within the quasiclasssical Green functions approach. A different 
approach which circumvent these difficulties while keeping a fully microscopic description of the problem can be 
obtained by expanding the field operators in a discrete basis and writing the Hamiltonian (1) in the form [Q 

H = - ^)4* c i<y + tiAoCjo + J2( A i 4i c h + A i c iT c a)> (2) 

where i,j run over the different sites of the system, tij are the hopping parameters connecting the different sites; \Xi 
and Aj being the chemical potential and order parameter in a site representation. The simplification introduced by 
this approach allows to deal with rather involved situations including spatial inhomogeneities (self-consistency) and 
non-stationary effects typically appearing in superconductors. For the voltage range eV ~ A the energy dependence 
of the transmission coefficients can be neglected and the transport properties can be expressed as a superposition of 
independent channels [ p"3|] . One can simplify even further the model to represent an SQPC with a single conduction 
channel, which can be described by the following Hamiltonian 

H = H L +H R + Y,^ {T) ' 2c l^R.a + t*e~^/2 c f Ra c La ) - fi L N L - m N R , (3) 

where Hl r are the BCS Hamiltonians for the left and right uncoupled electrodes characterized by a constant order 
parameters A^^ (for a symmetric contact A^ = Aji = A). 4>{r) is the time-dependent superconducting phase 
difference entering as a phase factor in the hopping terms describing electron transfer between the electrodes. In 
our model the transmission, a, can be varied between and 1 as a function of the coupling parameter t (see || for 
details). Within this model, the total current through the contact can be written as 

J ( T ) = | L ( <e#(T)/2 < c L(r)cB,(r) > ~ **e-^)/ 2 < ^ a (r)c L<7 (r) >) . (4) 

a 

The averaged quantities appearing in the current can be expressed in terms of non-equilibrium Green functions [ffBl . 
For the description of the superconducting state it is useful to introduce spinor field operators (Nambu representation) 
p6| , which in a site representation are defined as 



& = y ) - # = ( 4i «i ) . (5) 



-a 

Then, the different correlation functions appearing in the Keldysh formalism adopt the standard causal form 
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G% P (T a y f} )=-i<f$ i (T a )i>HT' f} )]> 



(6) 



where f is the chronological ordering operator along the closed time loop contour (jl5| . The labels a and f3 refer to 
the upper (a = +) and lower (a = — ) branches on this contour. The functions which can be associated within 

this formalism with the electronic non-equilibrium distribution functions [ p"7| , are given by the (2x2) matrix 



n_-f < cJ T (T / )c iT (T) > < Cj1 (t')c^(t) > 



<c] T (r')c l t i (r)> <c n (T>)cj l (T)> 



In terms of the G + the current is given by 



J(r) = -Tr 



a z (iG+- L (T,T)~PGi R (T,r) 



where t is the hopping element in the Nambu representation 

t 
1 -t* 



(7) 



(8) 



(9) 



The Green functions G~[~ are calculated using an infinite order perturbation theory with the coupling term in Eq. 
(3) considered as a perturbation. Within this approach these Green functions obey a set of integral Dyson equations 
[pj. As discussed in the next sections, the solution is strongly dependent on the biasing condition which determines 
the time dependence in the superconducting phase difference. 



III. CURRENT IN A VOLTAGE BIASED CONTACT 



The simplest biasing condition is that of a constant applied voltage. This situation is rather easy to achieve 
experimentally except for very small voltages (see section IV). In spite of its apparent simplicity the theoretical 
analysis is quite complex because of the time-dependent phase-difference which gives rise to a time dependent current 
containing all harmonics of the Josephson frequency u>q = 2eV/Ti, i.e J(r) = ^ n I n (V) expinwo- The current can be 
also decomposed into dissipative and nondissipative parts according to the different symmetry with respect to V of 
even and odd terms in the previous expansion 

In this case the integral Dyson equations can be transformed into a set of algebraic equations by a double Fourier 
transformation, defined by 

G„,mM = J dr J dT'e~ ZUJo{nT - mT "> /2 e tuj{T - T ">G(T,T') (10) 
An efficient algorithm for the numerical evaluation of the Green function Fourier components is discussed in Rcf. 

We shall concentrate in this section in the dc component of the current Iq which is the quantity more readily 
accessible experimentally. Fig. 2 shows the dc IV characteristics calculated from the fully quantum mechanical 
theory and from the semiclassical OBTK theory [Q . As can be observed the results become increasingly different for 
decreasing transmission. The fully quantum- mechanical calculation exhibits a pronounced subgap structure with steps 
at eV — 2 A / n which is hardly noticeable in the semiclassical theory. Both theories give the same result nevertheless 
for perfect transmission where interference effects, not included in the semiclassical theory, disappear due to the 
absence of backscattering. 

The experimental IV characteristics for atomic contacts of different metals are in remarkable agreement with our 
theoretical results. This is illustrated in Fig. 3 for the case of a one-atom contact made of Pb (these results are taken 
from Ref. (6|). This agreement has allowed to extract information on the conduction channels transmissions T n of 
metallic atomic contacts J5|,|6|,p^[ . 
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FIG. 2. dc current-voltage characteristics of a SQPC for different values of the normal transmission. Left panel corresponds 
to the semiclassical OTBK theory and right panel to the fully quantum mechanical calculation. 




voltage (A/e) 

FIG. 3. Measured current-voltage characteristics (symbols) for different realizations of a Pb one-atom contact at 1.5 K 
fabricated with the STM technique [6]. The full lines are numerical fits obtained by superposing four one-channel IV curves 
with different transmissions. 



The temperature dependence of the IV characteristics is shown in Fig. 4 for different values of the transmission. A 
remarkable feature of this dependence is that the SGS persists up to temperatures close to the critical temperature. 
When normalized to the temperature dependent superconducting gap the dc current exhibits a certain increase at 
low transmission, the opposite behavior being found close to prefect transmission. The crossover between these two 
tendencies is found for a ~ 0.8. 




eV/A(T) eV/A(T) 
FIG. 4. dc current-voltage characteristic for different temperatures and four values of the transmission. 



The limit of very small bias is particularly interesting due to the contribution of MAR processes of increasing order 
n ~ 2A/eV. This divergency in the V — > limit is eventually controlled by the presence of an inelastic relaxation rate 
77 (usually a small fraction of the gap parameter) which introduces a cut-off in the theory when eV < T). The effect of 
the inelastic relaxation rate is to damp MAR processes of order n > 2 A/77. As a consequence the system experiments 
a transition into a different regime where the total current becomes linear in V. In this regime the system response is 

determined by the adiabatic dynamics of the Andreev states at e(</>) = 1 — Tsin 2 ((/>/2) which move following the 
actual value of the superconducting phase. The total current can then be written as I(<j), V) = Is(<f>) + G((j>)V [ p0[ , 
where the supercurrent Is{4>) and the phase-dependent linear conductance G((j>) are given by 



G(4>) 



eA 



2h 



h I677 



yfl -a sin 2 (0/2) 
Aa sin 



tanh(^^) 



yjl -a sin 2 (0/2) 



sech(^M) 



(3V. 



(11) 



The expression of the supercurrent |14j2^^^] in Eq. (11) interpolates between the Josephson Is ~ sin0 behavior 
and the Kulik-Omelyanchuk Ig ~ sin 0/2 ballistic limit [ |23| . This behavior at high transmission has been recently 
confirmed experimentally using break junction techniques in an SQUID configuration jnj. It should be noticed that 
the expression of G(<j)) gives a definite answer to an old-standing problem concerning the form of this term known as 
the "cos0 problem" |24f| . The precise form of this term remains to be explored experimentally (a similar set up to 
that used in Ref. jy] could be used for this purpose). 

Finally, it should be stressed that from a mathematically point of view the two limits 77 — > and V — ► are not 
interchangeable || . In practice the limit eV — > with eV > rj can never be reached as there is always a finite although 
small inelastic relaxation rate present. 



IV. CURRENT BIASED CONTACT 



At very low voltages (and specially for high transmission) the contact impedance may become actually smaller 
than the voltage source impedance. If these conditions apply, the assumption of having an ideal source providing a 
constant voltage bias which fixes the phase dynamics is no longer valid. In this case one should take into account the 
electromagnetic environment of the contact in order to determine the phase dynamics and the system response to the 
external bias. 

For conventional tunnel junctions this limit is usually analyzed by means of the RSJ and RSCJ models |2p^| which 
represent the actual environment by a simple shunted circuit with a resistance R and a capacitance C connected in 
parallel to the junction. Within these simple models the phase dynamics is equivalent to that of a particle moving 
in a "tilted washboard" potential U((f>) = —It4> + / c cos0, where lb is the biasing current and I c is the Josephson 
critical current. At finite temperatures one should also consider thermal fluctuations acting as an stochastic force on 
the fictitious particle. 

To analyze the response of an SQPC under current bias one should generalize these models for contacts of arbitrary 
transmission |p5[. The description of the superconducting phase as classical variable will be valid as long as the 
Josephson coupling energy Ej ~ hl c /2e is much larger than the charging energy Eq ~ e 2 /2C. For an SQPC 
connected to a current source, the equations for the generalized RSCJ model would be given by 

V=^ (12) 

where is(0) and G(<f>) are given by Eq. (11) of the previous section and i n {(f>) is a fluctuating current whose power 
spectrum S is related to G(tf>) by the fluctuation-dissipation theorem S = AksTG (see section VI). It should be noticed 
that the above equations are strictly valid in the limit of small voltages induced on the contact, i.e. eV < 77, which is 
the condition for the validity of Eq. (11) in section III. The actual value of 77 is unknown but can be estimated to be 
of the order of A/100 or less pq| . In the mechanical analogy, the effective potential for the generalized RSCJ model 
can be written as 
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U(<f>) 



4e 
0h 



log 



cosh 



(13) 



and there appears a "position" dependent friction which comes from the dissipative term in G{4>). The inclusion 
of this phase-dependent term should have important consequences in the dynamics of the system. Notice that the 
particular form of G(4>) (Eq. 11) introduces a very asymmetrical friction with a minimum at the local minima of U((j)) 
and maximum at the local maxima. 

Integrating the Eqs. (12) of the generalized RSCJ model for arbitrary conditions is a formidable task. An ap- 
proximate solution for the overdamped case, i.e. G((f>)/C > (2eI c /HC) 1 ^ 2 , can be obtained following the procedure 
introduced by Ambegaokar and Halperin |27|| for overdamped tunnel junctions. The generalization of the Ambegaokar 
and Halperin theory is straightforward once we have identified the generalized potential (Eq. (13)) and the shunted 
resistance with 1/G((f>) (details will be given elsewhere). The measurement of the slope of the IV curve at zero 
voltage, which is directly related to G(</>), would provide information on the value of r\ in real systems. 



V. CONTACT UNDER MICROWAVE RADIATION 



As discussed in section III, the Andreev states play a central role in determining the adiabatic dynamics of an 
SQPC at low bias voltage. Considering that typical subgap energies are in the microwave range, it seems natural to 
propose using microwave radiation for a direct detection of Andreev states. This possibility has been suggested in a 
previous work by us (28) and in Ref. p9| ]. 

The effect of a microwave external field can be easily introduced in the single channel contact model. One can 
assume that the field intensity is maximum in the constriction region and neglect the effect of field penetrating inside 
the electrodes. Within this assumptions the field can be introduced as a phase factor modulating the hopping term t 
in Eq. (3), i.e. 

t{r) =te iaoCOSUrT , (14) 

where u> r is the microwave frequency, ao = eV op t / (huj r ) , V op t being the optical voltage induced by the field across the 
constriction. The parameter ao measures the strength of the coupling with the external field. The time-dependent 
hopping term can be expanded as 

t(r) =^z"J„(a )e m ^ T , (15) 

n 

where J n is the n-order Bessel function. For small coupling one can keep the lowest order terms in Eq. (15) and 
obtain some analytical results P^| . In the general case, the model Hamiltonian can be viewed, according to Eq. (15), 
as a superposition of processes where an arbitrary number of quanta of energy Tioj r are absorbed or emitted. As 
the temporal dependence of each term in Eq. (15) is formally equivalent to that in the constant voltage case, the 
generalization of the algorithm discussed in section III to the present case is straightforward. 
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FIG. 5. Induced dc current in a SQPC under microwave-radiation for different values of the transmission. The parameter 
oto controls the coupling to the external field (see text). 



Fig. 5 shows the induced dc current as a function of the microwave frequency for the case of low coupling constant 
(cvo = 0.1). All these results correspond to the situation in which the contact is carrying the maximum supercurrent. 
In this weak coupling limit the induced current is mainly due to the excitation from the lower to the upper Andreev 
state, which carries a negative current (i.e. opposite to the supercurrent). As a consequence the induced current 
exhibits a maximum for the resonant condition oj r = 2e (</>). At resonance, the induced current can be of the same 
order as the critical supercurrent. One can also notice a second stellite peak around e(4>) associated with two photon 
processes and a continuous band above A + e(<p). When the coupling constant cvo increases the contribution of higher 
order processes becomes progresively more important giving rise to a complex structure where the resonant condition 
for the excitation of the upper Andreev state can no longer be resolved p8[. 



VI. THERMAL AND SHOT NOISE 

The analysis of current fluctuations has a central role in the theory of transport in mesoscopic systems [ j30"t . 
Fluctuations can provide useful information on the microscopic dynamics (correlations) not contained in the average 
current. It is thus desirable to develop a fully quantum mechanical theory of current fluctuations in an SQPC on an 
equal footing as previously discussed for the current. In this respect, some limiting cases have already been analyzed 
in the recent literature: the excess noise for eV > A in a ballistic contact has been discussed in Ref. f3l|] , thermal 
noise for arbitrary transmission was analyzed in Ref. |j2[ while the case of perfect transmission and finite voltages 
has been addressed in Ref. J33|. 

The noise power spectrum is defined by 



S{oj, r) = fi dr' e lWT < SI(t + t')SI(t) + 5I{t)6I{t + t') > 



(16) 



where 81 (t) = J(t)— < I(t) >. For the evaluation of the above correlation functions a BCS mean field decoupling 
procedure can be performed. S(lu,t) can then be written in terms of nonequilibrium Green functions introduced 
in section III. In the voltage biased case, S(uj,t) can be expanded in harmonics of the Josephson frequency, i.e. 
S(uj,t) — ^2 S n (uj) exp inujQT. As in the case of the average current the noise Fourier components S n (u>) can be 
evaluated in terms of the Green functions matrix elements G n ^ m defined in section III. 

Let us start by analyzing the V = case where noise is due to thermal fluctuations. While in a normal QPC 
thermal noise has the usual well understood behavior, increasing linearly with temperature and with a flat frequency 
spectrum, in the superconducting case it exhibits very unusual behavior as a function of temperature, frequency and 
phase. Fig. 6 illustrates the frequency dependence of the thermal noise for different transmissions. As can be observed 
the noise exhibits two sharp resonances at oj = and oj — 2e(4>) corresponding to the excitation of the Andreev bound 
states. For ui > A + S exhibits a broad band arising from the continuous part of the single particle spectral 

density. 




cj/A 

FIG. 6. Power spectrum S(u>) of the zero-voltage current-fluctuations at UbT 
values. The contact is biased at the maximum supercurrent. 



0.2A in a SQPC for three transmission 
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The weight of the peaks at and 2e(0) can be evaluated analytically as discussed in Ref. [Q. We find that the 
zero frequency noise is related to the phase-dependent linear conductance by 5(0) = 4fcsTG ! (0) as expected from 
the fluctuation dissipation theorem. The exponential temperature dependence in G(<f>) gives rise to an exponential 
increase of thermal noise when fc^T ~ e(4>). It should be noticed that the ratio S(0)/2eIs(4>) can actually be divergent 
for any temperature provided that a — ► 1 and <fi — > n. On the other hand, the weight of the peak at 2e((f>) is zero 
for perfect transmission, increasing as a 2 (l — a), becoming the dominant feature for finite 1 — a and sufficiently low 
temperatures. In fact, the ratio between 5(2e(0)) and 5(0) is given by 



5(2e(0))/5(O) = 1(1 



\ 9 r 

a) tan —cosh 



<4>) 



k B T 



(17) 



Another quantity which is interesting to analyze and is directly amenable to experimental measurement is the 
shot-noise. Mathematically the shot-noise is given by the zero frequency dc component in the expansion of the noise 
power spectrum, i.e 5o(0), at eV ^> ksT. For simplicity we will consider the zero temperature case. Results for the 
shot-noise as a function of voltage are shown in Fig. 7 for several transmissions. The curves exhibit a pronounced 
subgap structure at the voltage values eV — 2A/n as in the dc current. In the case of the shot-noise the structure 
is more pronounced and is still observable for transmissions rather close to 1. In the perfect ballistic limit shot-noise 
is greatly reduced due to correlations associated with the Pauli principle as in the case of a normal ballistic contact 

0- 

On the other hand, in the tunnel limit the shot-noise is expected to reach the Poisson limit 5 ~ 2ql, where q is the 
transmitted charge in an elementary process. This relation offers the possibility to directly check whether multiple 
charges q — ne are actually being transmitted coherently in a n-th order MAR process |55|. Our theory allows to 
calculate the effective charges defined by the shot-noise current ratio. In the tunnel limit one finds that q exhibits a 
well defined step like behavior q/e = 1 + Int[2A/eV] confirming the above hypothesis p6[. 
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FIG. 7. Zero-frequency current-fluctuations at zero temperature (shot-noise) as a function of bias voltage. 



VII. SGS AND RESONANT TUNNELING 



The model discussed so far describes an SQPC with an energy independent transmission a. In some situations, 
that can be achieved experimentally, the normal transmission can have a non-negligible variation on an energy scale 
of the order of A. This can happen when the constriction region is weakly coupled to the electrodes by tunnel barriers 
as in the case of a small metallic particle or a quantum dot coupled to superconducting leads |$7|. We represent this 
situation by the following model Hamiltonian 

cr^Oa + cJ^Cw) + £o"-0<r + UfiQ^floi (18) 

where Hl and Hr describe the left and right leads, eo is a resonant level associated with the isolated constriction 
region, t v with v — L, R arc hopping parameters which connect the level to the left and right leads. The U term takes 
into account the Coulomb repulsion in the constriction region. The parameter U is basically the charging energy, Eq, 
and is related to the central region capacitance C, by U ~ e 2 /2C. For the subsequent discussion it is convenient to 
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introduce the normal elastic tunneling rates T v — 7r|tj,| 2 /3„(/x), where p v {p) are the normal spectral densities of the 
leads at the Fermi level. 

When the charging energy is much larger than both T and A Andreev reflections are completely suppressed and 
transport is only due to single-quasiparticle tunneling. This situation has been achieved in experiments on transport 
through nanometer metallic particles by Ralph et al. |57|. Model calculations presented by us in Ref. based on 
Hamiltonian (18) yield good agreement with the experimental results. 

We shall consider in more detail the case of small charging energy, in which the interplay between resonant tunneling 
and MAR give rise to novel effects and a very rich subgap structure p8p9| ]. Fig. 8 shows the dc IV characteristic 
for different positions of the resonant level eo with respect to the Fermi level. The tunneling rates are taken in this 
case as = = A. As can be observed, when the level is far from the gap region (case a) the limit of energy 
independent transmission is recovered and the subgap structure is similar to the one depicted in Fig. 2 (right panel). 
As the resonant level approaches the gap region the subgap structure becomes progressively distorted with respect 
to the energy independent transmission case. While the structure corresponding to the opening of odd-order MAR 
processes (i.e. at eV ~ 2A/n with odd n) is enhanced, the structure at eV ~ 2A/n with even n is suppressed. When 
r — » (not shown), one can also notice the appearance of resonant peaks in the IV characteristic for eV ~ 2eq. 




eV/A 

FIG. 8. dc current-voltage characteristic for a resonant level eo coupled to superconducting leads, eo = 5A (a), 2A (b), A 
(c), A/2 (d) and (e). Rn{0) is the normal resistance at the Fermi level. Curves (c), (d) and (e) have been displaced for the 
sake of clarity. 



VIII. CONCLUSIONS 



An overview of the results of a microscopic theory for the transport properties of superconducting quantum point 
contacts has been presented. These results include the response under different biasing conditions for an energy 
independent transmission as well as the case of resonant transmission. A remarkable agreement has been found 
between the calculated and the experimental dc TV curves for atomic-size contacts |^||. The agreement has allowed 
to extract information on the number and transmissions of the conduction channels in atomic-contacts of different 
metallic elements ||[l9] . On the other hand, the agreement shows the importance of interference effects included in a 
fully quantum-mechanical calculation and thus the need to go beyond semiclassical theories for describing this kind 
of systems. Additional predictions of the microscopic theory remain to be analyzed experimentally. For instance, we 
could point out the phase dependence of the linear conductance, the direct observation of Andreev levels in contacts 
under microwave radiation and the analysis of shot-noise. This last analysis would provide direct evidence of coherent 
transmission of multiple charges in MAR processes. Finally, the rich SGS in the presence of resonant transmission 
could be explored in S-2DEG-S devices wich are currently being developed (3). 
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